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Abstract 

The effect of a slowly rotating mass on a forced harmonic oscillator 
with two degrees of freedom is studied in the weak field approximation. It 
is found that according to the general theory of relativity there is a shift 
in the resonat frequency of the oscillator which depends on the density 
and rotational frequency of the gravitational source. The proposed shift 
is quite small under normal physical situations however it is estimated 
that for compact x-ray sources such as white dwarfs, pulsars, and neutron 
stars the shift is quite appreciable. 



1 Introduction 

The general theory of relativity is logically one of the most compelling theories of 
modern physical science. However being a physical theory it is imperative that 
it should not only describe as many different physical phenomenon as possible 
but also make accurate predictions of physical quantities relevant to observa- 
tions. As regards the accuracy in prediction the general theory of relativity 
has done extremely well in all its known tests such as the perihelion shift of 
mercury, deflection of light, and gravitational redshift experiments^'^. However 
the range of applicability where general theory of relativity can be compared 
with direct observations has remained small, owing to the fact that under nor- 
mal physical situations on earth and inside the solar system it is very difficult 
to differentiate its predictions from other rival theories such as the Newton's 
theory of gravitation, even when it is considered under its linear approxima- 
tion. The discovery of compact x-ray sources in the universe has opened up 

* Email address: bmmirza2002@yahoo.com 



1 



new avenues to the test general relativity in naturally existing systems. In these 
astrophysical laboratories extreme physical conditions are encountered which 
are impossible to be reproduced on earth. For example^ densities of a usual 
compact x-ray soince is of the order up to 10^* kg/m? and rotational periods 
IQ^^Hz. Under such conditions general relativity provides not only a consistent 
picture of various physical aspects of the phenomenon but also predicts, with 
unprecedented accuracy, magnitudes of various physical parameters involved 
such as the timing of rotational periods and their decay with time^. On the 
other hand all such systems are composed of atoms existing under a very strong 
gravitational field. Within the classical theory these atoms can be taken as 
forced harmonic oscillators whose physical behavior depends very much on the 
background gravitational field. 

In this paper we study the effects of gravitational field on the behavior of 
such oscillators within the framework of the general theory of relativity.. After 
giving a preliminary introduction, we write the equations of motion for a forced 
harmonic oscillator in section 2, making a plausible physical assiimption that 
the star is rotating slowly (i.e., far below the speed of light limit). In section 3 
we discuss the solution to the equations of motion and especially their behavior 
near the resonant frequency. In section 4, we give the behavior of the oscillators 
for compact astrophysical objects and also give a comparison. Lastly a brief 
summary and conclusion is given in section 5. 



2 The Equations of Motion 
2 . 1 Preliminaries 

According to the general theory of relativity^'^ the motion of a free particle of 
gravitational mass nig is given by the geodesic equation 

(fx°' ^„ dxl^ dx"' ^ ^ „ . „ o 
^+^^^^^ = «' -,^7 = 0,1,2,3 (1) 

where a:" = (x*,.t'^) = (r,t) are the position coordinates of the test particle 
and i = 1,2,3. Here s is the affine parameter and F^^ = g^"{gs0,'y + ff<57,/3 — 
5/37,5 )/2 is the Christoffel symbol, where g"'^ is the metric tensor and denotes 
differentiation with respect to x". 

Let the metric tensor be expressed as 

9oc0 = Vocp + hocH (2) 

where r;„^ = diag{—l, —1,-1,1) is the Minkowski metric tensor and hap is the 
perturbation to the metric. Then under the linear approximation ha/3 <C 1 we 
can take time t as the affine parameter. Requiring rotation to be slow (i.e., 
neglecting square and higher powers of velocity vector v =dv / dt) we obtain the 
following expression for acceleration for a slowly rotating sphere of homogeneous 
mass density 
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where 



(Pr 

^ = « + ^x« 



(3) 



G = - 



Vip, H = Vx4a 



(4) 



and 




'-dV, a = 




(5) 



p being the mass density and V is the vohime. Here the remarkable formal 
analogy of the above results with the classical electromagnetic theory must 
be noticed with the proviso that a factor of 4 is to be multiplied with the 
gravitomagnetic vector potential a. Pursuing this analogy, we find^ that for 
a slowly rotating sphere of homogeneous mass density first term in expression 
(3) is the gravitoelectric (GE) potential given by the Newtonian gravitational 
potential — Mf/r (where M is the mass of the gravitating source, and G = 1 = 
c unless mentioned otherwise). Moreover there is the gravitomagnetic (GM) 
potential v x H where 



R being the radius and Q, is the angular velocity of the gravitational source. 
The second part of equation (3) has a non-Newtonian origin. It is regarded 
to exhibit typically Einsteinian effects and therefore plays an important role 
in testing Einstein's theory of gravitation in the weak field and slow rotation 
approximation^'^. The independence of the GM potential from a particular 
frame and from a particular coordinate system used has been demonstrated 
thus making it physically significant^*^. This effect can be interpreted as 'grav- 
itomagnetic current' induced in the vicinity of the gravitational source due to 
its rotation. For its analogy with the Lorentz force law for the electromagnetic 
field the force m,g{G + v x H) is often referred to as gravitoelectromagnetic 
(GEM) force. However to measure the effects of this force on a given physical 
system high accuracy in experiments is required where it is necessary to isolate 
the experimental setup from contingencies^^'^^. 

2.2 Forced Harmonic Oscillations in Weak Gravitational 
Field and Slow Rotation Approximation 

Let us consider a forced damped harmonic oscillator having two degrees of 
freedom placed in a GEM field. If mj is the inertial mass of the test particle 
then in addition to the binding force — m/Wgr , the damping force — m/7V, and 
an external force Fq exp iu)t acting along, say, the x direction, there also acts on 
the particle a GEM force mg(G+v x H). Considering the GE field of magnitude 
—g to be along y-axis and the GM field of magnitude H in the direction normal 




(6) 
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to the xy- plane we find that the GEM force has a component rUgHvy in the 
X direction and a component —mgHvx in the y direction; where Vx and Vy are 
the components of velocity in the x and y directions respectively. With these 
additional terms and the fact that m.g = m.i{= m), in accordance with the 
general theory of relativity, we get the equations of motion as a coupled system 
of ordinary differential equations 

(Px 2 dx Fo . ^^dy 

— ^ = -ujqX + — e^^iujt + H—, (7) 

dt'^ dt m dt 

d'^y 2 dy dx 

Now let us introduce the auxiliary variable ^ = x + iy, then the equations 
of motion give a single equation 



—^ + {^ + iH)—+wl'St=—e-^^iwt-ig (9 
dt-^ dt m 

Notice that the GM parameter H has a coupling to the damping parameter 

7 in the complex plane {x,y). Dimensional analysis of equations (7) and (8) 

shows that H must have dimensions of cycles per second i.e., of frequency. 

These considerations give rise to the question whether these 'gravitomagnctic 

oscillations' produce significant effects on the physical behavior of an oscillatory 

system. As we shall see in the following that the answer to this question must 

be given in affirmative. However these effects are extremely small except in 

those cases where density and rotational frequency of the gravitational source 

are extremely large, such as those found in compact astrophysical x-ray sources. 



3 Solution to the Equations of Motion and In- 
terpretation 

The exact solution to equation (9) can be obtained, however we are interested in 
studying the behavior of amplitude especially near the resonant frequency. To 
do so we express 5R as A exp itot where A is the complex amplitude. Substitution 
in equation (9) gives 

— — in cxp(—iLut) 

M^) = — — — 10 

Multiplying the last expression with its complex conjugate and taking square 
root we have the following real expression for the amplitude 



I,, / (fo/m)2 + g2 _ 2(fo/m)ff sina;i 

Comparing this expression for the amplitude with the expression for the 
amplitude of a forced damped harmonic oscillator having two degrees of freedom 
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we notice that there are two main differences, one the addition of g dependent 
terms in the numerator, two an addition of the term —ujH, in the denominator, 
to the square of eigenfrequency wq of the system. The second difference is of 
particular importance since it involves the effect of GM field on the resonant 
frequency of the system. To estimate this effect we notice that in equation (11) 
for small damping i.e., 7 <C w and given values of the parameters Fo/m ,g, 
and ujQ, the amplitude A is maximimi when a;§ — — uiH ~ i.e., when the 
applied frequency is approximately {—H±^yH^~i^4uj^)/2. Expanding in powers 
of {H/2oJo)'^ and requiring that H <^ uq we obtain the following approximate 
result for resonant frequency 

a;~wo-Y, 7<1, H <§: ujq (12) 

This result shows that there is a shift in the resonant frequency by an amount 
H/2. 



4 Application to Compact Stars 

The outer atmosphere of a white dwarf consists mainly of fully ionized *i?e, ^^C, 

and ^^O atoms forming a crystalline structure in which electrons can be taken as 
forced harmonic oscillators, moving under the influence of lattice vibrations^'*. 
Along the radial direction the motion of these electrons is very much limited 
by the GE attraction and the Fermi pressure, therefore in any thin layer en- 
veloping the star the electrons can be taken as oscillating with two degrees of 
freedom. As the star of mass greater than the Chandrasekhar limit exhausts 
its thermonuclear fuel the GE force not only contracts the core (increasing the 
mass density per unit area), its angular frequency also is increased from IHz 
(= 1 cycle per second) to IMHz or above. To study the resonant behavior of 
these electrons we begin with relatively less severe case of white drawf of mass 
IMg, = 1.989 X 10^"% and radius 7 x lO^m rotating at frequency IHz . Then 
the case of the pulsar of mass 1.4M® and radius 3 x 10*m rotating at frequency 
30Hz is studied. Finally we consider a neutron star of mass 2M(gand radius 
1 x 10**TO rotating at frequency IkHz. In these cases g (= GM/R^ , where 
G = 6.67 x IQ-'^'^Nm? /kg"^) ranges from 2.7075 x lifm/sec^ (for the white 
dwarf) to 2.0636 x lO^m/sec^ (for the pulsar) to 2.6533 x W^m/ sec^ (for the 
neutron star); consequently {Fo/m)'^ + 9^ + 2 {Fq /m)g sin ivt varies from order 
lO^^m/ sec^ to lO^^m/sec^ where 6 = ojt is the phase. Now from equation (6) 
we have in the equatorial plane at the surface of the star i. e., r = R : 

H^\H\=^^^\n\ (13) 

where fj, = 4,G/bc^ = 0.5940 x 10~^'^m/kg . For the above given values of the 
relevant physical parameters we have from expression (13), H is 0.1687x 10~^i?2; 

, 1.6540i7z, and 236.2932_ffz for the whitedwarf, the pulsar and the neutron 
star respectively. With these values of parameters we plot the resonance curves. 
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using expression (11), for the three cases in figure (1) and (2) for electron with 
eigenfrequency IkHz ^ 9 = 9i = 0, and 7 ~ O.lHz. 

We notice that for a slow rotating pulsar the resonance shift Aw =| lu — luo \ = 
H/2 is not very large (about 0.8269iJ2) as compared with the shift for the case 
of the neutron star (about 117.3196i?z). 

5 Summary and Conclusion 

The above considerations can be summed up as follows: under the influence of 
gravitomagnetic field of a rotating compact star there is a shift in the resonant 
frequency of an oscillator lying in its vicinity. The shift depends on the density 
and on the rotational frequency of the star. 

Inside the solar system and especially in terresrial laboratories the detection 
of gravitomagnetic resonance shift is very unlikely (for Earth H is approximately 
0.6440 X lO^^^Hz and for the sun it is 0.6597 x 10~^^Hz). However for compact 
gravitational objects the shift lies in the domain of observations. These obser- 
vations may not be directly possible at present, but with further theoretical 
as well as observational developments, it may well become possible to provide 
evidence for such an effect. 
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5.1 Figure Captions 

Figure 1: Amplitude plots for an oscillating test particle exhibiting the resonance 
shift for a white dwarf (mass 1M@ , radius 7 x lO^m , rotational frequency 
liJz)and a pulsar (mass 1.4Afg), radius 3x lO^m , rotational frequency 30Hz) as 
the gravitational source. The frequency oj is measured in Hz and the amplitude 
is in units m/Hz. sec^. 

Figure 2: Amplitude plots for an oscillating test particle exhibiting the res- 
onance shift for a pulsar (mass 1.4M(S), radius 3 x lO^m , rotational frequency 
30Hz) and a neutron star (mass 2Mg, , radius lO^m , rotational frequency 
IkHz) as the gravitational source. The frequency uj is measured in Hz and the 
amplitude is in units m/Hz. sec^. 



7 
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